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Cross section analyses in MiniBooNE and SciBooNE 

experiments 

Teppei Katori for the MiniBooNE and SciBooNE collaborations 

Massachusetts Institute of Technology, Cambridge, MA 

Abstract. The MiniBooNE experiment (2002-2012) and the SciBooNE experiment (2007-2008) are modem high statistics 
neutrino experiments, and they developed many new ideas in neutrino cross section analyses. In this note, I discuss selected 
topics of these analyses. 
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THE MINIBOONE EXPERIMENT 

The MiniBooNE experiment uses the Booster Neutrino Beamline (BNB) [1] and the MiniBooNE detector [2]. The 
BNB creates < Ey >~700 (600) MeV muon neutrino (anti-neutrino) beams. The neutrinos (anti-neutrinos) travel 
~520 m before observed by the MiniBooNE Cherenkov detector filled with mineral oil (CH2), through productions of 
charged particles. 
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Signal definition 

How to define a signal event is important for any cross section measurement. We encountered a problem when the 
charged-current 1 7r+ production (CCl;r+) to charged-current quasielastic (CCQE) cross section ratio was studied [3]. 
The final state interactions (FSIs) affect the pion observation, through pion absorption, charge exchange, pion pro- 
duction, and re- scattering. All of them have large errors. How to remove these FSI effects to measure genuine pion 
kinematics from the neutrino interaction vertices? The answer we came was not to remove FSI effects, but define our 
signals differently. Since the FSIs are not understood well, any corrections on FSIs introduce extra biases in the data 
and make the data model dependent. Instead, we define the signals from the final state particles in the detector. In the 
case of the CCl;r+ interaction, a signal event is defined by "1/i and l7r+ in the final state". 

Figure 1 shows this situation. Fig. la shows a CC process accompanied with one pion production, where both the 
muon and the pion leave the nuclear target and are observed successfully. 

However, some final-state pions may also be created not by the primary neutrino interaction, but by FSI processes, 
such as the hadronic re-interaction (Fig. lb). Since our detector cannot distinguish such a pion from those pions created 



a) Signal 






d) Signal 



pton 

absorption 

in detector media 




FIGURE 1. Signal definition of pion production measurement, (a), (b), and (d) are classified to be "signal" 
Notice both (c) and (d) are pion absorptions, but (b) is in the target nuclei, and (d) is in the detector media. 



and (c) is "not signal" 
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FIGURE 2. MiniBooNE CCItt+ZCCQE total cross section ratio [3]. In the left plot, there is no correction apphed about FSI 
in the target nuclei, and hence, reactions are called "CCl;r+-like" and "CCQE-like". In the right plot, the simulation dependent 
correction is applied. 

from the primary interaction, such event must also be included in what we define as "signal". The data is a sum of 
both, which help theorists to study FSIs. 

Figure Ic shows an opposite case, where pions made by the primary neutrino interaction fails to leave the target 
nuclei, hence are not observed. Traditionally, experiments estimate how many pions disappeared due to FSI, and apply 
corrections. However, such corrections are model dependent and should be avoided. By our signal definition, this kind 
of event is classified as "not signal". 

The last example (Fig. Id) is a cumbersome situation. Pions can be absorbed during their propagation in the detector 
medium. Since this is a detector dependent effect, we need to correct based on our simulation. This introduces an 
error, and for example, such detector dependent nuclear error dominates the MiniBooNE CC l7r° production (CCl;r°) 
measurement [4] . 

The signal definition affects how the final data appear. Fig. 2 shows the CCl;r+/CCQE total cross section ratio 
taken from Ref. [3]. These 2 total cross section ratio plots are based on different signal definitions. On the left, signal 
is defined from observables (Fig. 1), instead of primary interactions, and hence the measured total cross section ratio 
is called "CCl;r-like" to "CCQE-like" ratio. This result is independent of any nuclear models. The price to pay is, 
theorists can compare their models with our data only if they have ways to simulate FSI effects. Otherwise, theorists 
are recommended to use the right plot. Here, simulation dependent corrections of the target nuclei FSI are applied, 
therefore the ratio can be interpreted as a primary interaction cross section ratio. Historically only such measurements 
are performed [5,6]. The price to pay is, the resulting ratio is nuclear model dependent, and the error bars are inflated 
to take into account this model dependency. 



Data driven correction 

Measured interactions always include background events. And these background events need to be removed, 
however, this depends on the predictions from the simulations. To avoid such model dependency, MiniBooNE actively 
uses the data driven corrections to the background events, except for the CCl;r+ cross section measurement where the 
signal purity is ~90% [7]. 



CCQE cross section measurement 



In the CCQE cross section measurement [8], CCl;r+ events make the biggest background. If a CCItt^ event loses 
a pion through FSI, its final state particles are identical with those of CCQE events. Simultaneous measurements of 
the CCQE and CCl;r+ candidate sample are performed, and information from the CCl tt^ candidate sample allows 
the background in the CCQE candidate sample to be corrected as a function of reconstructed 4-momentum transfer 
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FIGURE 3. Examples of MiniBooNE data driven corrections. In the left plot, dirt event enhanced NCEL sample is shown with 
various components from the simulation [9]. The template fit corrects the scaling factor of dirt events. In the middle plot, the 
ratio of measured CCl7r+ to simulated CCItt^ double differential cross section is shown [4]. This is used to correct the CCl;r+ 
background for CCl;r° cross section measurement. In the right plot, flux scale factors of neutrino component in the anti-neutrino 
beam are shown as a function of Ey [13, 14]. There are 3 independent measurements and they are consistently lower than 1 (default 
simulation). 



(2^). This corrected background is subtracted from the CCQE sample to measure various cross sections. Note this 
correction is valid only within the precision of FSI. To allow theorists to study FSIs, MiniBooNE also published the 
subtracted background distributions [8]. 



Neutral current elastic (NCEL) cross section measurement 

In the NCEL cross section measurement [9], backgrounds created outside of the detector (dirt events) are significant. 
Fortunately, the MiniBooNE detector is big enough to see the spatial profile of dirt events in the detector. Figure 3, 
left, shows the dirt event enhanced NCEL sample as a function of the vertex location in Z (beam direction) [9] . The 
template fits can find the normalization factor of dirt events, as a function of measured total nucleon kinetic energy. 
After scaling the dirt events, data and simulation agree well. This technique to correct dirt events is also used for the 
oscillation analysis [11]. 



CCl n"" cross section measurement 

In the CCl7r° cross section measurement [4], CCl;r+ events make large backgrounds. Unlike the CCQE cross 
section measurement where the CCItt^ background is corrected as a function of reconstructed Q^, here the correction 
is based on the double differential cross section of pion kinetic energy and neutrino energy (Ey). Figure 3, middle, 
shows the data to simulation ratio [4], where data is the measured CCl;r+ double differential cross section [7]. This 
ratio is applied as a correction to the simulated CCl;r+ background. 



Neutral current 1 7r° (NCItt"^) cross section measurement 



In the case of the NCl7r° cross section measurement [10], no data driven correction is applied. However, the 
measured rate of NCl7r° is used to correct the background distribution of the oscillation candidate events [11]. Here, 
7r° background in the oscillation sample is corrected as a function of 7r° momentum [12]. This measurement not only 
corrects the rate, but also constrains the 7r° rate error, because the cross section error of NCI ;r° is ~25% from the 
simulation, whereas the measured NCl7r° rate has ~5% error. 
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FIGURE 4. The MiniBooNE CCQE signal sample based on different background removing methods [19]. The red lines are for 
background subtraction method, and the blue lines are for signal fraction method. 



Anti muon-neutrino CCQE cross section measurement 

The V^CCQE cross section measurement [13] has a complicated structure, because of the presence of neutrino 
contamination in the anti-neutrino beam. This is measured, and corrected by using 3 independent techniques [13, 14]. 
Fig. 3, right, shows the result. 

The CCl7t~ production from anti-neutrino interactions make large backgrounds for the v^CCQE cross section 
measurement. The simultaneous fit technique used in the V^CCQE measurement cannot be applied since the CC\7t~ 
final state contains same final state particles with v^CCQE (7r~ is ~100% nuclear captured). Therefore the CCl;r~ 
distribution is corrected by assuming the same kinematic distribution as the CCln^ distribution in neutrino mode. 

Finally the v^CCQE cross section model is corrected based on the relativistic Fermi gas model tuning [15] from 
the v^CCQE measurement [8]. Note that the absolute cross section measurement does not depend on the cross section 
model of the signal channel (in this case, v^CCQE model), except in unfolding. 

The v^NCEL cross section measurement [16] is underway. With this result, IVLiniBooNE will complete all 4 
quasielastic and elastic scattering cross section measurements. The difficulty is, our NCEL cross section measurements 
are not on a proton or a neutron target, but the combination of both. This makes it difficult to apply interesting ideas 
to access to nucleon parameters [17, 18]. Theorists are encouraged to invent new ways to utilize ]V[iniBooNE CCQE 
and NCEL data! 



Background removing process 

Background events need to be removed from the data, but depending on the knowledge of the estimated background, 
there are several ways to remove the background. 

If the knowledge of the background events includes the absolute scale (for example, backgrounds are measured in 
situ), the background subtraction method is applied {dt — hi). In this case, signal events in /th bin is simply a difference 
of /th bin of data minus /th bin of background. This is a simple and preferred method, because it does not depend on 
the simulation of the signal events. The drawback is, you may get negative bins. 

If the knowledge of the background is at most a fraction of the total event, the signal fraction method is applied 
(4 X ^q^)- Ii^ this case, signal events in /th bin is data times predicted signal events (^0 divided by predicted total 
events {st + bt) of the simulation. The drawback is, the signal sample becomes signal model dependent. 

Figure 4 shows the example of these two methods from the IVLiniBooNE v^CCQE analysis [19]. Since the data- 
simulation agreement is good in entire kinematic region, the difference of these two methods is small in this sample. 



Unfolding 

Unfolding is an important stage of the analysis, because measured kinematics are often smeared or distorted 
by detector effects. These detector effects need to be corrected by the data unfolding. IVLost of IVLiniBooNE cross 
section data are unfolded by the iterative Bayesian method [20], to avoid the fast oscillation problem which is often 
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FIGURE 5. The 7r° kinematics unfolding from NCl7r° cross section measurement [22]. Left plot shows n° momentum, and the 
right plots is for 71° angular distribution. 

seen for histograms with many bins unfolded by the inverse matrix method [21]. The unfolding technique depends 
on many factors, and every single histogram needs to be assessed for the best unfolding scheme. Figure 5 shows 
NCl7r° kinematic distributions based on different unfolding techniques [22]. In this analysis. Four different unfolding 
techniques (one of four is no unfolding) are used to compare results. In the left, Tikhonov regularization [23] is chosen 
to unfold 71° momentum, but the same technique does not work for 71° angular distribution (right), and the iterative 
Bayesian method is chosen for the main result. 

THE SCIBOONE EXPERIMENT 

The SciBooNE detector is an X-Y tracker at the BNB [1]. The detector consists of 3 parts, the Scintillation bar vertex 
detector (SciBar, consists of CgHg), electron catcher (EC), and muon range detector (IMRD) [24]. Inclusive CC and 
CCQE total cross section results show similar excesses to the IMiniBooNE cross section results [25, 26]. 



Event classification 



SciBooNE can classify each event based on the topology. Figure 6 left, shows the CCItt^ analysis flow chart of 
event classification [27]. An event is classified depending on the type of tracks, number of tracks, and amount of vertex 
activity (energy deposit around the neutrino interaction vertex). Fig. 6 right plot shows the event sample with "1 muon 
and 1 pion tracks" in the final state, where the lower vertex activity region has more coherent CCl TT^ production 
according to the simulation. However, the data is consistent with the absence of coherent pions, as first observed by 
the K2K experiment [28]. The vertex activity is a powerful parameter, and similarly, NCl;r° analysis confirmed the 
presence of coherent production by the same technique [29] . This information is used for the latest T2K oscillation 
analysis [30]. 

The high resolution of the SciBar detector allows detailed study of each recorded track. The azimuthal angular distri- 
bution of pions may reveal further information of coherent pion production mechanism [31]. Single proton momentum 
measurement from the NCEL scattering [32] allows neutrino energy reconstruction without lepton kinematics. SciBar 
can also tag complicated topologies, such as CCl;r° [33]. All these are possible due to the high resolution of the SciBar 
detector. When more tracks are measured in more detail, an event can be classified in even smaller sub divisions. This 
eventually allows one to study the detailed structure of the FSIs. In this conference, ArgoNeuT showed an event-by- 
event counting of protons from CC interaction. Future high resolution experiments, such as IMicroBooNE [34], will 
allow further classification to understand more details of neutrino interactions. 



SclBar-MRD matched sample 



w/ activity 



PRD-5toppe<| 



Number of 
tracks 




MftD-penetrated 



Somii selection 



• DATA 

I I CC coherent n 

^§ CC resonant n 

■ Other 

■ CCQE 



^3thTvt^ w^octivity w/ooZ;]^ \/ 




20 30 

Energy deposit (MeV) 



FIGURE 6. Flow chart of event classification for the CCl 7r+ analysis (left). MRD stopped, 2 tracks, 1 muon-like and 1 pion-like 
sample is shown on the right plot. 
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